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Abstract 

We derive a finite set of nonlinear integral equations for describing the 
finite size dependence of the ground state energy of the 0(4) nonlinear sigma 
model. By modifying the kernel functions of these equations we propose 
nonlinear integral equations for the finite size effects in the SS-model. The 
equations are formulated in terms of two complex valued unknown functions 
and they are valid for arbitrary real values of the couplings. 



1 Introduction 



The SS-model is a two-parameter family of two-dimensional integrable quantum field 
theories (QFT), which possesses U{1) x U (1) symmetry and its action can be written 
in terms of three boson fields with an exponential interaction [1 . In this theory the 
scattering matrix is the direct product of two S-matrices of the sine-Gordon (SG) 
model with different coupling constants. This family of models includes as particular 
cases some other interesting QFTs like the 0(4) nonlinear sigma (NLS) model, the 
anisotropic principal chiral field and the N = 2 supersymmetric SG model. 
The action of the model can be expressed by three boson fields 0: 

Ass = J fx l^jydM'' + ^ [cos(a0i + a02)e^'^« + cos(a0i - 602)6"^'^^] | , 

(1) 

where the parameters a, a and (3 satisfy the relation : 

+ - /^^ = 47r. (2) 

The S-matrix of the theory is : 

s{e) = -s,{e)®Sp{e), (3) 

where Sp{9) is the SG S-matrix with |^ = The relations between p,p and the 
parameters of action are : 

Action is unitary if > so there are restrictions on the parameters p, p : 

p + p>2. (5) 

If j9 > 1 and p > 1 the spectrum of the model consist of four particles with equal 
masses. If p ^ oo and p ^ oo then the S-matrix Q becomes SU{2) x SU{2) 
invariant and it will coincide with the S-matrix of the 0(4) NLS model. Thus the 
SS-model can be regarded as a two-parameter f/(l) x U{1) invariant deformation of 
the 0(4) NLS model. 
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Once the spectrum and the associated S-matrix is known, the TBA method can 
be employed for calculating the ground state energy of the model in finite volume. 
This was achieved j]P for the SS-model, and for integer values of p and p the TBA 
equations can be encoded in an extended (afiine) "Pp+p Dynkin-diagram (see figure 
E). For general values of p and p the form of the integral equations are much 
more complicated, and the number of unknown functions depends on the continued 
fraction form of p and p j,2.^. In the asymptotically free case {p,p —>■ oo) the number 
of unknown functions becomes infinite in the TBA equations. The solutions of the 
TBA equations also satisfy the well-known Y-system equations jlj. 




Figure 1: Dynkin-diagram associated with the SS-model Y-system 

In p] it was noted that if a model has an S-matrix in the form of a direct product 
Sg ® Sh, and the TBA equations are already known for the models described by 
S-matrices So and Sh, and the individual TBA equations are encoded on Dynkin- 
like diagrams of type G and H respectively, each having one massive node, then the 
TBA equations for the model with the direct product S-matrix can be obtained by 
gluing the individual TBA equations together at the massive node. This method 
can also be applied for the SS-model due to the tensor product form of the S-matrix. 

Another approach for calculating finite size effects in a QFT is the use of some 
integrable lattice regularization of the QFT which can be solved by the Bethe Ansatz 
method. In models that can be solved by Bethe Ansatz, the well-known T- and Y- 
systems 121111 naturally appear and from them, using the analytical properties of the 
T-functions, one can easily derive the same TBA equations which can be obtained 
from the S-matrix. Using such an integrable regularization, an alternative method 
was worked out by Destri and de Vega [7j for calculating finite size effects in the 
SG model. This is called the nonlinear integral equation (NLIE) technique. The 
advantage of this method is that one gets a single integral equation for any real 
value of the coupling constant, and the method can be extended to excited states, 
too jH El ^1- The problem is that this method can only be applied when the 
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ground state of the model is formed by real Bethe Ansatz roots. In most cases this 
is not true. For example the ground state of the spin-S" XXX-chain is formed by 
25-strings O IH 113 CHI • 

Recently J. Suzuki JI] managed to derive some new types of nonlinear integral 
equations for describing the thermodynamics of the higher spin XXX-model at 
finite temperature, which can be regarded as the generalization of the spin-^ case 
which was treated successfully earlier (ZlEDj- These new equations can be regarded 
as a particular mixture of TBA and NLIE, which gives back the p ^ oo (isotropic) 
limit of the Destri-de Vega equation [7j in the spin-^ case. 

Last year motivated by the gluing idea |3] and Suzuki's results [TT] Dunning pro- 
posed nonhnear integral equations [l2j similar to Suzuki's equations for describing 
the finite-volume ground state energy of some perturbed conformal field theories 
whose S-matrices can be written in the form of a direct product. These equations 
can also be regarded as a mixture of TBA and NLIE, but these equations have the 
advantage that the number of unknown functions is less than in the TBA equations, 
hence these equations are much more convenient for numerical studies. Another 
advantage is that if the model has a single coupling constant then the equations are 
valid for all real value of the coupling constant. 

In this paper we derive a finite set of nonlinear integral equations for describing 
the finite size dependence of the ground state energy of the asymptotically free 0(4) 
nonlinear sigma model. The equations are formulated in terms of two complex valued 
unknown functions and hence they are more convenient for numerical calculations 
than the infinite set of TBA equations ^11^. By modifying the kernel functions of 
these equations we propose nonlinear integral equations for the finite size effects in 
the SS-model for finite values of the couplings. The equations are also formulated 
in terms of two complex valued unknown functions and they are valid for arbitrary 
real values of the couplings. The form of these equations can be regarded as an 
appropriate gluing of two SG NLIEs together. 

The paper is organized as follows. In section 2, we recall the main results of the 
light-cone lattice approach to the 0(4) NLS model. In section 3, we introduce T- 
and Y-systems from the Bethe ansatz solution of the higher spin six-vertex model. 
In section 4, using the light-cone lattice approach and Suzuki's [TT] method we derive 
nonlinear integral equations for describing the finite-volume ground state energy of 
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the current-current perturbation of the SU{2)2s WZW model. In section 5, using 
the hght-cone lattice approach and motivated by Suzuki's method we derive new 
nonlinear integral equations for the same model. In section 6, in the S* — > oo limit 
we glue together the two nonlinear integral equations to a much simpler equation. 
This equation describes the finite- volume ground state energy of the 0(4) NLS 
model. In section 7, we propose nonlinear integral equations for describing the finite- 
volume ground state energy of the SS-model by modifying the kernel functions of 
the equations obtained for the 0(4) NLS model. In section 8, we perform some 
analytical and numerical tests on the proposed equations. In section 9, some special 
cases are considered. The summary and conclusions of this paper are given in section 



2 The light-cone lattice approach to the 0(4) NLS 



Our starting point is the Bethe Ansatz solution of the integrable lattice regular- 
ization of the 0(4) NLS model [221 123] • In this section we briefly summarize the 
results of this approach j2Sj • The fields of the regularized theory are defined at sites 
( "events" ) of a light-cone lattice and the dynamics of the system is defined by transla- 
tions in the left and right light-cone directions. These are given by transfer matrices 
of the isotropic higher spin six- vertex model with alternating inhomogeneities. This 
approach is particularly useful for calculating the finite size dependence of physical 
quantities. We take N points (A^ is even) in the spatial direction and use periodic 
boundary conditions. The lattice spacing is related to /, the (dimensionful) size of 
the system : 

21 

(«) 

The physical states of the system are characterized by the set of Bethe roots 
{xj, j = 1, . . . , M}, which satisfy the Bethe Ansatz equations (BAE) 



10. 
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where 

M 

k = 2S, Q{x) = Y[{x-Xj), (8) 

i=i 

and Xq is the inhomogeneity parameter. 

The energy (E) and momentum (P) of the physical states can be obtained from 
the eigenvalues of the light-cone transfer matrices: 

if (g+P) ^ / 1 xAl Qi^O ~ '^^) i^{H-^P) ^ / -. nm Q(~^0 + /q\ 

^ ^ Q(a;o + ^A;)' ^ Q^-^^ ^ ^k)■ 

Besides the usual procedure, taking the thermodynamic limit (A^ — > oo) first, fol- 
lowed by the continuum limit (a 0) one can also study continuum limit in finite 
volume by taking N ^ oo and tuning the inhomogeneity parameter xq simultane- 
ously as 

2 4 2 2N 

xo = - log = - log — , (10) 

vr ma ir ml 

where the mass parameter m is the infinite volume mass gap of the theory. If we take 
this continuum limit at a finite, fixed spin value S we get the energy and momentum 
eigenvalues of the current-current perturbation of the SU{2)2s WZW model. After 
this continuum limit we shall take the S ^ oo limit in order to obtain the energy 
and momentum eigenvalues of the 0(4) NLS model in finite volume with periodic 
boundary conditions |221- 

It is important to realize that this procedure does not give all the eigenvalues of 
the Hamilton and the momentum operators of the 0(4) NLS model! Representing 
the 0(4) symmetry of the model as 0(4) ~ SU{2)l x SU{2)ji, one can show that 
only those states of the Hilbert space of the theory can be described by this method 
which are SU{2)l (or equivalently SU{2)r) singlets [23 El- Since the ground state 
of the 0(4) NLS model lies in this sector of the Hilbert space, this description is 
appropriate to examine its finite size dependence. We note that the same results 
could be obtained for the 0(4) NLS model by using the approach of refs. [T^ IT3]. 



3 Fusion hierarchy and T-systems 

In this section we introduce T-systems from the fusion hierarchy of the isotropic 
spin-S* six-vertex model. Let Vi ^ C''^^ be the irreducible SU{2) representation 
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with spin /j/2, and let Rfj'^^\x) be a linear operator acting on Vi (S> Vj. These 
operators depend on a complex parameter x and satisfy the Yang-Baxter equation 

R?^'''\x)R^^''\x + y)Rtt'\y) = Rt''\y)R^lt\^ + y)R^u''\x). (n) 

The solutions of this equation can be obtained by fusion ^7] from the simplest 
R-matrix which corresponds to the spin-i representation both in Vi and V2. This 
R-matrix is of the form 

3 

=a; + z + «^a''®(T", (12) 

a=l 

where a" are the Pauli matrices. The concrete form of the R'^^^'^^\x) matrices can be 
taken from ref. ^Hj. From these R- matrices one can define the family of monodromy 
matrices with alternating inhomogeneities Xi = (— l)*"'"^Xo: 

T^P\x, {x,}) = \x -xi-t{k+p- 1)) . ..R^^^'\x -XN-t{k + p- 1)), (13) 

where k = 2S. These matrices act on the tensor product Va Vh where Vh = 
Vi ® ... ® Vtv is the quantum space of the system and Va is the auxiliary space with 
~ C'^, \4 ~ C^. In the rest of the paper we consider the case when is even, 
because this is necessary for the light-cone lattice approach. Define the transfer 
matrix by taking the trace of the monodromy matrix over the auxiliary space; 

Tp{x,{x,}) = TTaT^P\x,{Xi}). (14) 

Due to the Yang-Baxter relation (fTT|) the transfer matrices ()14|) form a commutative 
family of operators acting on Vh IT^ 

[rp(x, {Xi}), Tniy, {Xi})] = 0. (15) 

Due to the integrability and commutativity all these transfer matrices can be si- 
multaneously diagonalized by the algebraic Bethe Ansatz jTH]- The eigenvalues of 
the transfer matrices can be characterized by the solutions of the Bethe Ansatz 
equations [THl 



{xj + Xo — ik) {xj — xq — ik) j Q{xj — 2i) 



N 



{xj+xo + ik){xj-xo + ik)\2 _ Q{xj + 2i) ^ ^^^^ 
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The eigenvalues of the transfer matrices p4|) are of the form 

Q{x + i{p+l)) Q{x + \ 
Q{x + i{2l - J9 - 1)) Q{x + i{2l - p - 3)) ' 



T(r\-Sr w/^'Vrl Q{x + l{p+l)) Q{x-l{p+l)) 

lp^X)-2^Pl W [Li) 



1=1 
where 

p-l 1-2 

= Yl (p{x + 2im-i{k+p-l))Yl4>{x + 2in + i{k + l-p)), (18) 

m=l—l n=0 

N N 

= (x — Xo)^(x + Xo)^. (19) 

The eigenvalues ()17j) of the transfer matrices (fT^ satisfy the so-called T-system 
equations: 

Tp{x + i)Tp{x -i) = fp{x) + Tp_i(a;)Tp+i(x), (20) 

where 

p 

fp{x) = n '^(^ + '(2J -p + k)) 4>{x + i{2j -p-k-2)). (21) 
i=i 

From a T-system ()2U|) one can define a Y-system as follows 

= ^'-''^'^'K (22) 

y,W = l+,,W = S(^±i§i^. (23) 
These functions satisfy the Y-system equations jH IH] 

yj{x + i)yj{x + i) = Yj.i{x)Yj+i{x). (24) 
One can redefine the T-system elements with 

Tj {x) fj (x) = aj {x)Tj (x) , (25) 
where (Tj{x) satisfies the relation 

aj{x + i)<Jj{x — i) = (Tj_i(x)crj+i(x), (26) 
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then the new Tj{x) functions satisfy the same T-system relations as (pUj) but with 
different fj{x) functions: 

fji^) fji^) = (^ji^ + i)(^jix - i)fj{x). (27) 

This is called gauge transformation. Under such a transformation the Y-system (j221- 
is invariant. As we will see later the energy of the model can be expressed by an 
element of the gauge invariant Y-system. On the other hand one has the freedom 
of choosing that gauge for the T-system (j^IH) which is convenient for the particular 
calculations. 



4 Nonlinear integral equations I. 

In this section we recall the derivation of Suzuki's equations jlT] to describe the 
ground state energy of our model in finite volume. Consider the following gauge 
transformation of the T-system pn|l 

fp{x + i)fp{x -i) = fp{x) + fp_i(,T)fp+i(x), (28) 



where 



and 



U^) = \{<p{x + i{p-k- 2j)) <P{x -i{p-k- 2j)), (29) 



T_i(a;)=0, To{x) = l. (30) 

The solutions of these equations are of the form |TTj 



where 



p+i 

Tp(x) = 5^A!^)(x), (31) 



^ ^ ~ ^' Q{x + i{2l -p-l)) Q{x + i{2l -p-3)y ^^'^ 



'^^P\x) = Yl (p{x + i{p-k- 2j + 1)) JJ</)(x -i{p-k- 2j + 1)). (33) 
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We define the following auxiliary functions JT] 

*W = S-i^fM, , = 1....,, (34) 

F,(x) = l + y,(x) = ^H^±il^H^, J = l,...,fc (35) 

+ = l + (37) 

From ()28H33|) it follows that the auxiliary functions ()34ll37p satisfy the following 
functional relations: 

n{x + z)=p_<P{x + 2zj)^^^^B{x), (38) 

Ux-z) = l[<j>{x-2zj)^^±^B{x), (39) 
j=i 



(x) Q{x -ik ~ 2i) ^ 

nL0(a;-22j) ~g(^+W 



Ka:) = ^.^^^ .... (41) 



fi(x)fi(x) = Y,{x), (42) 

yj(x + i)?/j(x-i) = Fj_i(x)Fj+i(x) j = l,...,A;-2, (43) 

yk-i{x + i)yk-i{x + i) = Yk-2{x)B{x)B{x), (44) 

fk-i{x + i)fk-i{x -i) = fk-i{x)Yk-i{x). (45) 
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We introduce two other auxiliary functions 

^i{x) =Q{x-ik), <i!2{x) =Qix + ik). (46) 

In order to be able to derive integral equations from these functional relations one 
needs to know the positions of the zeroes and the poles of the auxiliary functions 
()34II37|) . Due to the relations ()38H45|) we only need to know the zeroes of Q{x) 
and Tj(x). For the ground state of our model the zeroes of Q{x) form y pieces 
of /c-strings, but for finite values there are deviations from the string hypothesis 
|24j . From [21] one can see that those Bethe roots have the smallest deviations from 
the imaginary value of the string hypothesis that have imaginary parts ±{k — 1) 
according to the string hypothesis. These deviations are always less than 1/2 and 
this is important because the analytic properties of ^1^1(0;) {'^2{x)) are influenced 
mainly by these roots on the upper (lower) half plane of the complex plane near the 
real axis. 

The transfer matrices Tj{x) j = l,...,k ()28|) have no zeroes in the ground 
state in the "main" strip < |Imx| < 1. We can list the strips where the auxiliary 
functions are analytic and non zero (ANZ). 





^i{x) 


ANZ 


Im a; < 1/2, 












^2(a;) 


ANZ 


Im X > -1/2, 










b{x] 


,Bix) 


ANZ 


< |Im x\ < 1/2, 










b{x] 


,B{x) 


ANZ 


< |Im x\ < 1/2, 














ANZ 


< Im x < 1, 


J 


= 1,. 


.,k- 


1, 




Yjix) 


ANZ 


< Im x\ < e. 


j 


= 1,.. 


.,k- 


1 






ANZ 


< Im X < 1, 


J 


= 1,- 


.,k. 





We introduce new variables by shifting the arguments of 6(x),i3(x) and h[x),B[x) 
by ±^7 mi 

ao(x) = h{x — i-y), Uo{x) = B{x — 27) = 1 + ao{x), (48) 

ao{x) =b{x + i'y), Uo{x) = B{x + i'y) = 1 + ao{x), (49) 
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where 0<7<l/2isan arbitrary real parameter. This shift is necessary because the 
original functions b{x),b{x) have zeroes on the real axis. Due to the ANZ property 
of Tk{x) ()47p and the fact that lim|^|_>oo ^ ^ogT/^^x) = the following relation holds 

m 

oo oo 

= [ dx4- logTfcfx - z) e'«("-') - ! dx4- logTfcfx + i) e'«("+'). (50) 
J ax J ax 

— oo — oo 

From ()38|) and ()39|) one can express Tk{x ± i) with the auxiliary functions and by 
substituting these expressions into (j^Uj) one gets in Fourier space (see conventions 
for Fourier transformation in appendix A.) 

7, ^ / N ,x i.„ smhikq) cos(xoa) 

dl^r{q>Q) = mNe-^'^ 1;{ . 

cosh(gj smh^gj 

e(l-7)9 £-(1-7)9 _ 

2 cosh(g) 2 cosh(g) 

dWiiq < 0) = 0, (52) 

dW2{q > 0) = 0, (53) 

J/v,,(,<0) = -mNe'^ ""^[^^^ 

cosh(g) smh(g) 

£(1-7)9 _ _ g-(l-7)i? _ 

2 cosh(g) 2 cosh(g) 

where we introduced the notation 

oo 

dlF{q) = I dx e^^^ ^logF(x). (55) 



One can derive similar relations for dlyj{q)'s and dlYj{q)'s from ()43|) . and dlTk_i{q) 
and dlYk-i{q) from Substituting these relations into the definitions of ao(a;) 

and ao(x), one obtains /c + 1 algebraic relations in Fourier space. After taking 
the inverse Fourier transformation of these relations and integrating over x, we get 
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Suzuki's hybrid nonlinear integral equations ^1] 



log?/i(x) = {K *\ogY2){x), 

\ogy,{x) = (K*logF,_i) + (ir*logF,+i)(x), j = 2, . . . , A; - 2, 

\ogyk-i{x) = {K*\ogYk.2){x) + {K+^*\ogUo){x) + {K-^*\ogUo){x), 

logao(x) = P^(x-27) + (F*log?7o)(x)-(F+2(i->)*log[/o)(x) (56) 

+ (ir-^*logn_i)(a;), 

logao(x) = V^{x + t^) + {F*logUo){x)-iF-''^'-^U\ogUo){x) 

+ {K+^*logYk.,){x), 

where {K* f){x) = J dy K{x — y) f{y) is the convolution and the "source" function 
T>n(x) on the lattice reads as: 



T^n{x) = iA^arctan 



sinh 



■K{x+i) \ 
2 ) 



cosh (^) 

The kernel functions of ()56|1 are of the form 

K{x) 



2, f2N\ , , 



4 cosh(7ra;/2) ■ 



(5^ 



dq 



\q\~iqx 



27r 2 cosh(g) 



and we have used the notation 



(59) 



f^^{x) = f{x±tri). 



(60) 



We get the continuum limit of eqs. by taking the — oo limit. The ker- 

nel functions do not change because they are independent of N, but the "source" 
function changes in the continuum as 



/ 7TX \ 

I){x) = lim V^ix) = — m/cosh ( — ) 

N^oo V 2 / 



(61) 



In the TBA language the complex auxiliary functions aQ{x) and ao{x) resum the 
contributions of those Y-system elements whose index is larger than k — 1. Equations 
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fl56|) are graphically represented in figure El The big "bubble" denotes the complex 
auxiliary functions which resum the contributions of those TBA nodes, which are 
inside it. In our notation the names of the complex unknown functions and the 
kernel function are indicated. 



o — o o — p O 

1 2 k-1 \kk+lk+2 

ao{x), ao{x), F(x) 

Figure 2: Graphical representation of nonlinear integral equations I. 



The energy and momentum of the model can be easily expressed by our 
auxiliary functions due to the fact that apart from some trivial normalization factors: 



e^K^^+P) ^ T,(xo + 0, and e^t^^"^) ~ n{-x, - z). 

After some straightforward calculations [3 ESj one gets in the continuum limit for 
the ground state energy of the model 

oo 

' Tl{x — 27) 



Eq{1) = Ebuik ~ -J I cosh 



\ogUo{x) 



(62) 



^ fdx cosh (^i^) log t/o(.) 



where Ehuu. is an overall divergent factor in the continuum and it has the form on 
the lattice: 



E, 



bulk 



-^^^ 1 



21 ^ I 

j=0 



2(1 + 2j) - 0:0 
i{l + 2j) + xo 



when k is even, (63) 



iV2 

Ebulk — I X 



(fc-3)/2 



(2xo)+ 7 



log 



j=0 



22(1 + j) -Xq 



when k is odd, 



where 



X[x) = / — 



dq sin(ga;) e 



-kl 



27r g 2 cosh(g) 



(64) 
(65) 
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From these formulas one can easily see that Eh^ik is independent of the Bethe roots, 
thus independent of the dynamics of the model and so Ebuik appears in the energy 
expression as a bulk energy term same for all states of the model. One can easily see 
that Ehuik is divergent in the continuum, but after subtracting this divergent term 
from the energy the finite integral terms in ()(i2|l describe the finite size dependence 
of the ground state energy of the model. 

5 Nonlinear integral equations II. 

In this section we derive a new set of nonlinear integral equations to describe the 
ground state energy of our model in finite volume. Although these equations are 
similar to Suzuki's equations they differ in that the "source" term is coupled to the 
TBA part of the system and in the TEA language the complex auxiliary functions 
resum the contributions of those Y-system elements whose index is larger than k. 
(See figure El) Consider another gauge transformation of T-system 




(66) 



where 




(67) 





The solutions of these equations are of the form 




(69) 



1=1 



where 




(x) = To{x + i{2l - p - 2)) 



Q{x + i{p + 1)) Q{x - i{p + 1)) 



(70) 



Q{x + i{2l -p-l)) Q{x + i{2l - p - 3)) • 



We define the following auxiliary functions: 





(71) 



14 



r,(x) = l + y,(x) = ^H^±^^H^, j = l,...,A: + l, (72) 

Mx) = '^''^^"''^,^;;'^'^^"~'\ n{x) = 1 + Mx). (74) 

From ()fifill7()j) it follows that the auxiliary functions ()71II74|1 satisfy the following 
functional relations: 

f,+i(x + z) = fo(x + i{k + 2)) ~ ~ ^'j 7^(x), (75) 

r,+i(x - = ro(x - ^(fc + 2)) + + n{x), (76) 

C^(x — ik — i) 

/i(x) = -— — Tfc(x), (77) 

^ ^ To(x + z(A; + 2)) g(x - ^A; - ^) ^' ^ ^ 

^ ^ To(x-^(A; + 2)) Qix + ik + t) ^' ^ ^ 

7^(x)7i(x) =n+i(x), (79) 

yj{x + i)yj{x - i) = Yj_i{x)Yj+i{x) j = 1, . . . , - 1, (80) 

2/fc(a; + i)yk{x + i) = Yk^i{x)n{x)n{x) , (81) 



fkix + i)fk{x -i) = foix + i{k + 1)) fo(x - i{k + l))n(x). (82) 

In order to be able to derive integral equations from these functional relations one 
needs to know the positions of the zeroes and the poles of the auxiliary functions 
()71II74|) . Due to the relations ()75ll82p we only need the zeroes of Q{x) and Tj{x). The 
zeroes of Q{x) are the same as in the previous section, but the Tj(x)s have changed 
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due to the gauge transformation. This change is only an overall factor which does 
not depend on the Bethe roots. So in this gauge all Tj{x) can have y-fold degenerate 
zeroes in the "main" strip at places ±xo, ixo ± i, but these zeroes cancel from the 
Y-system elements fl71H72|) except the yk{x) case which has y-fold degenerate zeroes 
in the "main" strip at ±xo. These zeroes will give the standard TBA source term in 
our final equations 123] . Now we can list the strips, where the auxiliary functions 
are analytic and non zero (ANZ): 

Im X < 1/2, 
Im X > -1/2, 
-3/2 < Im X < 0, 

< Im X < 3/2, (83) 
< |Im x| < 1, j = l,...,k-l, 
< |Im x| < e, j = j, . . . , k e > 0, 
< |Im x| < 1. 

We introduce new variables by shifting the arguments of /i(x),7Y(x) and h{x),'H{x) 
by ±27': 

a(x) = /i(x-i7'), U{x) =n{x -i'y') = l + a{x), (84) 



^i(x) 


ANZ 


^2(x) 


ANZ 


h{x),T-C{x) 


ANZ 


h{x),'H{x) 


ANZ 




ANZ 




ANZ 




ANZ 



a{x) =n{x + 1^), U{x) =n{x + = l + a{x), (85) 

where < 7' < 1/2 is an arbitrary real and fixed parameter. This shift is necessary 
because the original functions 7i(x),7i(x) have zeroes and poles on the upper and 
lower half plane respectively. Due to the ANZ property of T^+i (x) (j83|l and the fact 
that lim|2;|_»oo ^ logTfc+i(x) = the following relation holds 

00 00 

0= [ dx -^logffc+i(x-i) e^''^^-^) - [ dx -^logffc+i(x + 2) e^''(^+^). (86) 
J dx J dx 
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From (fTSj) and (fTBj) one can express Tk+i {x ± i) with the auxihary functions and by 
substituting these expressions into (jH?)|) one gets in Fourier space 

dm{q>0) = ^ !y! + :--r^ dlU{q) - - dlU{q), (87) 

2 cosh(g) 2 cosh(g) 2 cosh(g) 

dWi{q<0) = 0, (88) 
dW2{q > 0) = 0, (89) 

c^/vl>2(g<0) = - . /;V - T^-rn ^^^(g) + ^^-[TTT ^^^^g)' (^0) 

2cosh(g) 2cosh(g) 2cosh(g) 

where 

dlf{q) = dlf-^^^^\q) - e"' dlf^^^^^\q). (91) 

We do not need the exphcit form of dlTQ^''~^'^\q) for our calculations, but we only 
have to know that the following identities hold 

e" dlf^'^^^^\q) = dlf~'-''^'\q), e'" dlf^^''^^\q) = dlf^^''^^\q). (92) 

After a similar procedure that has been done in the previous section one gets the 
following nonlinear integral equations in the continuum 

log vi (x) = {K* log Y2) (x) , 

\ogy,{x) = (K*logF,_i) + (ir*logr,+i)(x), j=2,...,A;-l 

\ogyk{x) = V{x) + {K*\ogYk-i){x) + {K+^' *\ogU){x) 

+ {K-^' *\ogU){x), (93) 

loga(x) = (F*logf/)(a;)-(F+2(^-^')*logf/)(x) + (fs:-^'*logFfc)(x), 

loga(x) = (F*logf/)(x) - (F-2(^-^'Ulogf/)(x) + (is:+^' *logFfc)(x), 

where the "source" function and the kernel functions are the same as in the previous 
section ()58ll6ip . In the TBA language the complex auxiliary functions a{x) and d{x) 
resum the contributions of those Y-system elements whose index is larger than k. 
(See figure El) 

The energy of the continuum model can be expressed by a gauge invariant Y- 
system element: 

00 

Tfi f / TTX \ 

Eo{l) = Ehuik -- dxcosh (^— j logn,(x). (94) 
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a{x), a{x), F{x) 

Figure 3: Graphical representation of nonlinear integral equations II. 



Equations (j56|) and (j93|) are different descriptions of the ground state energy of 
the current-current perturbation of the SU{2)]. WZW model. In this paper one 
of our main purposes is to describe the finite size dependence of the ground state 
energy of the 0(4) NLS model with a closed finite set of nonlinear integral equations. 
According to the light-cone lattice approach the A; — > oo limit of equations (j^Uj) 
and (jnni) describes the ground state energy of the 0(4) NLS model. In this limit 
the equations ()56p and ()93p contain an infinite number of unknown functions, but 
as we will see in the next section we are able to construct a finite set of nonlinear 
integral equations by combining ()5f)|l and ()93|1 in the k oo limit. 



6 Nonlinear integral equations for the 0(4) NLS 
model 

In this section we construct a finite set of nonlinear integral equations from 
and ()93|1 in the k —>■ oo limit for the ground state of the 0(4) NLS model. The 
complex auxiliary functions 1)40114111 and ()73II74|1 are gauge dependent and have a 
finite continuum limit only in an appropriate gauge, but the Y-system elements 
()22II23)1 are gauge invariant and have finite continuum limit. One can see from (jl^ 
and (|7^ that some gauge invariant Y-system elements can be expressed by these 
non-gauge invariant auxiliary functions: 

t/o(x + z7)f/o(x-z7) =n(x), (95) 



U{x + i-f')U{x - i-f') = Ffc+i(x). (96) 

If we choose the Y-system elements (1221123)1 (j = 1, 2, . . . ) as auxiliary functions, we 
get the standard TBA equations of the model, which can be encoded in an infinite 
diagram (figure EK-). 
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O O O • O O a O O • O O b 

1 2 k-l k k+1 k+2 k-2 k-1 k k+1 k+2 

Figure 4: a.) Dynkin-diagram associated with the current-current perturbation of 
the SU{2)k WZW model, b.) Dynkin-diagram associated with the 0(4) NLS model. 

In the A; — > cxD limit the diagram will be symmetric with respect to the massive 
node (figure EbO; so in the A; — > oo limit Yfc_i(x) = Yfc+i(x). So in the A; — > oo limit 
the following relations will also be true 

Uo{x + ij)Uo{x-i-f) =Yk{x), (97) 



U{x + i-i')U{x - i-i') = Ffc-i(x). (98) 

Substituting (j98|) into the equations for ao{x),ao{x) in (j56|l and (j95|l into the equa- 
tions for a{x),a{x) in ()93|) and deforming the integration contour appropriately we 
get a finite set of nonlinear integral equations in the k oo limit, which corresponds 
to the 0(4) NLS model case. The equations are as follows: 

log ao{x) = V{x - Z7) + {F* log Uo) (x) - {F+'^^^-^'^ * log Uo) (x) 

+ (K+(^'"^) * logf/)(x) + (ir-(^'+^) * \ogU){x), 

logao(a;) = P(a; + «7) + (F*logf/o)(x)-(F-2(i-T)*log?7o)(x) 

+ (K+(^'+^) * log U) (x) + (ir+(^-^') * log U) (x) , 

loga(x) = (F*logf/)(a;) - (F+2(^-^') *logf7)(x) (99) 

+ (K+(^-^') * logf/o)(x) + (K-(V+7) ^ logf7o)(x), 

loga(x) = (F*log[/)(x) - (F~^(^~^'Ulogf/)(x) 

+ (ir+(^+^') * logf/o)(a;) + (ir+(^'-^) * \ogUo){x), 
Uo(yX) = l + ao{x), Uo{x) = 1 + ao{x) , U{x) = l + a{x), U{x) = l + d{x). 

Equations (jU^ are graphically represented in figure 

The energy expression is the same as (jHSI)- This is a closed set of equations for 
four complex unknown functions "a(x), a(x), ao(x), ao(x)" , but a{x) is the complex 
conjugate of a{x) and aQ{x) is the complex conjugate of ao{x), therefore only four 
real unknown functions describe the model. 
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o — o a o — o — q — •- 

fc-3 k-2 k-iy \k k+1 k+2 J \ k-3 k-2 k-Xj k V^'+l '^'+2 



a(x) , a(x) ^ F ix) ao{x), ao(x), F{x) a(x), a(x), F{x) a{x) , a{x) , F (x) 

Figure 5: Graphical representation of the nonhnear integral equations of the 
0(4) NLS model 

From one can derive a different set of nonlinear integral equations to describe 
the finite size dependence of the ground state energy of the 0(4) NLS model. Using 
fl98|) in and deforming the integration contour appropriately we get the following 
equations 

\ogyo{x) = V{x) + 2{K+"'' *\ogU){x) +2{K~"'' *\ogU){x), 
loga(x) = (F*logf/)(x) - (F+2(i-V) ^logf7)(a,) 

+ {K-^' *\ogYo){x), (100) 

loga(a;) = {F *\ogU){x) ~ {F-^''^-"''U\ogU){x) 

+ (ir+^'*logFo)(a;), 

U{x) = l + a{x), U (x) = 1 + a{x) , 

where we introduced the notations: 

yo{x) = lim yk{x), Yo{x) = lim Yk{x). (101) 

k~*oo fc— >oo 

Equations ()1()()|1 are graphically represented in figure Eb- In this case the ground 
state energy of the 0(4) NLS model is of the form 

oo 

Eo{l) = Ehuik - J / t^a;cosh (^— j logyo(a:). (102) 

— oo 

Equations ()100|1 contain only three real unknown functions because yo{x) and Yo{x) 
are real and a{x) is the complex conjugate of a{x). 

Although equations pOOj] are more convenient for numerical calculations in the 
0(4) NLS model than equations (jU^ but, as it will be seen in the next section, it is 
a generalization of that describes the finite size dependence of the ground state 
energy of the SS-model for finite couplings. 



20 



7 Nonlinear integral equations for the SS-model 



In this section we propose nonlinear integral equations to describe the finite size 
dependence of the ground state energy of the SS-model. In equations the 
kernel function F{x) apart from some trivial factor, is nothing but the derivative of 
the logarithm of the soliton-soliton scattering amplitude at the = Svr point, 



l~/7rT\ ~ 

F{x) = -Goo (— ) , G^e) = -^logsue) 



(103) 



The scattering matrix of the SS-model is a direct product of two SG S-matrices 
with different couplings. The 0(4) NLS model is a special case of the SS-model and 
its scattering matrix is a direct product of two SG S-matrices at the jS"^ —>■ Stt point. 
In equations (jU^ we recognize a building block of the S-matrix of the 0(4) NLS 
model, as the kernel function F{x). As we mentioned in the introduction the SS- 
model can be regarded as a two-parameter deformation of the 0(4) NLS model and 
we assume that the nonlinear integral equations, which describe the finite size effects 
in the SS-model, can also be obtained by a two-parameter deformation of equations 
()99|). The deformation is given by changing the kernel function F{x) Gp{x) 
with appropriate p value, where Gp{x) is basically the logarithmic derivative of the 
soliton-soliton scattering amplitude: 

oo 

G{x)= [ ^ e-^^^ sinh((p-l)g) ^^^^^ 
^ J 2n 2 cosh(g) sinh(pg) ' 



Gp{x) = \g, (^) , GM = -^-^logsue) 



(105) 



"p+1 
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V{x - Z7) + (Gp * \ogU,){x) - (G+2(i-7) * logf/o)(x) 



The conjectured equations are of the form: 
log ao(a;) 

+ (ir+(^'-^) * \ogU){x) + {K-^'<'+'<) * logf/)(x), 

= P(x + «7) + (Gp * logf/o)(x) - (G;2(i-^) * logf/o)(x 

+ (ir+(^'+^) * log f/) (x) + (ir+(^-^') * log U) (x) , 

= (Gp_i * log U) (x) - (Gt_^^"'^ * log [/) (x) 

+ (/^+(^-^') * logf/o)(x) + (/^-(^'+^) * logf/o)(x), 



log ao(x) 
loga(x) 
loga(x) 



(106) 



(Gp_i * logf/)(x) - (G.2f-^') * \ogU){x) 



'p-i 

+ (ir+(^+^') * logf/o)(x) + (ir+(^'-^) * logf/o)(x), 
= l + ao(x), f7o(x) = 1 + ao(x), f/(x) = 1 + a(x), ?7(x) = 1 + a(x), 

where p and p are the parameters of action (1) and 0<7<l/2, 0<7'<l/2 
arbitrary, fixed real parameters. Equations ()106|) are graphically represented in 
figure EK- 





Ti '^^"jSx^-i 



a(x),a(x),Gp_i(x) aQ(3;),aQ(x), Gp {x) a(3;),a(x), Gp„ i (x) a(3;),a(x), Gp _ i (x) 

Figure 6: Graphical representation of the nonlinear integral equations of the 
SS-model 

The energy expression is 

oo 

— CX) 

(107) 

Due to the relations a(x) = a*(x), ao(x) = a*(x) (* denotes complex conjugation) 
we have only four real unknown functions for all real values of the two couplings 
p,p. The form of equations (jlOfij) is valid for p > 1 and p > 2, because in this case 
the kernel functions of ()106|) have no poles in the "main" strip (0 < Imx < 1). 

In contrast to equations (jnTljl and (jlOOj) . which were based on a derivation, equa- 
tions ()106|1 are only conjectures. For this reason in the next section we will give 
some analytical and numerical evidence showing that equations ()106p are indeed 
describing the finite size effects of the ground state energy in the SS-model. 
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8 The test of the equations 



In this section we will make some analytical and numerical tests on our conjectured 
equations ()l()fi|) . First we calculate the UV central charge of the model using the 
equations ()106|) . Using the standard method of refs. [201 1211 the energy in the 
conformal limit (/ 0) can be expressed by the dilogarithm functions. The energy 
in the conformal limit is of the form: 

Eo{l) ^ (^L+(a+(oo)) + L+(a+(oo)) + L+(ao+(oo)) + L+(ao+(oo)) (108) 
-L+(a+(-oo)) - L+(a+(-oo)) - L+(ao+(-oo)) - L+(ao+(-oo))^ , 
where L^{z) is defined by the integral 

z 

M--) = -y<(.(^L_J-— j, m.) = l(_), (109) 



and L{x) is Roger's dilogarithm function. The functions a^{x), a+(x), ao+(x), ao+{x) 
denote the kink functions corresponding to a{x), a(x), ao{x), ao{x) respectively. The 
limits of these kink functions at infinity follows: 

a+(oo) = a_|.(oo) = 1, a+{—oo) = a+(— oo) = 0, (HO) 



ao+(oo) = ao+(oo) = 0, ao+(-oo) = ao+(-oo) +oo. (Ill) 
From these using the simple identities: 



L+(oo) = L(l) = -, L41) = -, (112) 
one gets the ground state energy in the / — >• limit 

^o(0 ^ -f . (113) 

From this one can easily see that the effective UV central charge of the model is 
equal to three {cuv = 3). We obtained this value from our conjectured equations 
fll06p and it agrees with standard TBA calculations Ij. So we have checked our 
equations pO(j|l analytically in the small / limit in leading order. 
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Next we will check our equations in the large / limit analytically in next to 
leading order. The equations (jlOfij) can be solved iteratively in the large / regime. 
After some easy calculations one gets that the ground state energy is of the form 

Eq{1) = E(i) + ^(2) + ©(e-^^'"'), ^^'^ = Ef^ + Ef\ (114) 

where 

oo 

^(1) = I dd cosh(0) e-'"'=°^^(^), (115) 

2vr J 



— oo 
oo 



= ^ I de cosh(^) e-2-'-^h(e)^ ^^^g) 



47r 



= y rf^ cosh(e) e-'"'=°^'^('') y d0' $(0 - ^') e- 

— oo — oo 

where the kernel function ^{6) is of the form 



■ml cosh(6') 



(117) 



<^{9) = / did <l(cj) e-^"^ $(cj) = <lp(cj) + l'p(cj), (118) 



where 

sinh2(^) sinhf^^4^ 

$ (cj) =4-8 . — ^ ^ (119) 

' cosh(^) sinh(£f^) ^ ' 

The same results can be obtained from the TBA equations of the SS-model with 
integer parameters p,p using the method of [2ZI. One can recognize that ^{9) is 
nothing but 1/i times the trace of the logarithmic derivative of the two-body S- 
matrix of the SS-model: 

^(^) = -^-a Tr2log5(e), S{e) = -{S, ® S^m. (120) 
I do 

It is well-known that the ground state energy in finite volume is related to the free 
energy density of the system at temperature T by /(T) = Eo{l/T)T. The free 
energy density can be expressed as a series in the low temperature limit (virial 
expansion) and the coefficients of the series are completely determined by S-matrix 
data ^H] • If we transform the next to leading order large / expansion of the ground 
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state energy ()1 14111 19j) into the next to leading order low temperature expansion of 
the free energy, we get perfect agreement with the prediction of the virial expansion 

m- 

So far we have made some analytical tests on our proposed equations ()106|) in 
the small / and large / regimes. In the intermediate regime we solved numerically 
the equations ()106|) for some integer values of p and p, and we have solved the 
corresponding TBA equations, too. In every case we examined perfect agreement 
between the two sets of numerical data was found. 



9 Special cases and forms 

In this section we consider first the special case p = p. In this case instead of ()l()(i|l 
it is better to use the following modified version of (jlOOj) 

log?/o(x) = V{x) + 2{K+''' *\ogU){x) +2{K-^' *\ogU){x), 

loga(x) = (Gp_i*logf/)(x)-(G+5'~^'Ulogf/)(a;) 

+ {K-^' *\ogYo){x), (121) 

loga(x) = (Gp_i*logf/)(x)-(G;5'-^%logt/)(x) + (ir+^'*logFo)(x), 

U{x) = l+a{x), tJ{x) = l + a{x) Yo{x) = 1 + yo{x) . 

The graphical notation of eqs. 1)1211) can be seen in figure Eb- The energy is of the 
form 




— oo 



These equations are valid for the p > 2 case and they are better for numerical 
calculations than ()106|) because the number of unknown functions is less than in 

(Unni). 

Next we consider the special form of equations ()in6)l . which corresponds to the 
N = 2 supersymmetric SG model. This model is a special case of the SS-model. 
The scattering matrix of the model is of the form [SH] 

= -iS2^S,m, (123) 

/32=87rp 
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where is the couphng constant of the SUSY-SG model. We can see that the S- 
matrix of this theory is a special case of the S-matrix (jS} of the SS-model. Applying 
equations ()106|) for this case and taking 7' — we obtain 

logao(a;) = P(a; - ^7) + (Gp * logf/o)(x) - (G+^^i-^) * logt/o)(x) 

+ {K-^ *\ogU){x) + {K-^ *\ogU){x), 

logao(a:) = P(a; + ^7) + (Gp * log f/o)(a;) - (G^^d-^) * fog f/o)(a;) 

+ (fs:+^*logf/)(x) + (ir+^*log[7)(x), (124) 

loga(a;) = {K^^ *\ogUo){x) + {R-^ *\ogUo){x), 

\oga{x) = {K+^ *\ogUo){x) + {K-"' *\ogUo){x). 

Uo{x) = l + ao{x), f7o(a:;) = 1 + ao(a;), U{x) = l + a{x), tJ (x) = 1 + a{x) . 

It can be seen that in this case a{x) = a{x) G M, and we get the same equations as 
the ones proposed by Dunning |T21 for the N = 2 SUSY-SG model. 



10 Summary and Conclusions 

From the isotropic higher spin six-vertex model we derived two different sets of non- 
linear integral equations describing the finite size dependence of the ground state 
energy of the current-current perturbation of the SU{2)k WZW model. In the 
k ^ 00 limit we glued the two set of equations together and obtained a finite set 
of nonlinear integral equations for the 0(4) NLS model. Starting from these equa- 
tions we proposed a new set of nonlinear integral equations for describing finite size 
effects of the SS-model. We made analytical and numerical tests on our conjectured 
equations and all the results of these tests made us confident that our conjectured 
equations are correct. The advantage of our equations is that the number of un- 
known functions is minimal, and the equations can be defined for all real values of 
the coupling constants. 

It would be interesting to generalize these equations for describing (all) excited 
states energies of the model. 
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A Appendix 

In this appendix we summarize our conventions of Fourier transformation. Let 
f{x) be a bounded continuous function defined on the whole real axis. Its Fourier 
transform is defined as: 

oo 

f(q) ^ j dx e"^^ fix) (125) 

— oo 

and the inverse Fourier transform is 



oo 




(126) 



— oo 



The convolution of two functions is defined as follows: 

oo 

{f*9){x)^ J dy f{x-y)g{y). (127) 

— oo 

The Fourier transform of such a convolution is the product of the Fourier transforms 
of the two functions. 
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